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Abstract 

On the basis of conditions for existence of nuclear burning wave in neutron- multiplicating 
media (U-Pu and Th-U cycles) we show the possibility of surmounting so-called dpa- 
parameter problem, or in other words, an algorithm of optimal nuclear burning wave 
mode is proposed yielding parameters of the wave (fluence/neutron flux, width and speed 
of nuclear burning wave) satisfying the dpa-condition, characterizing the tolerable level 
of radioactive stability of reactor materials, in particular the cladding materials. 

It is shown for the first time that capture and fission cross-sections of ^^^U and ^^^Pu 
increase with temperature within 1000-3000 K range, that under certain conditions may 
lead to global loss of stability of nuclear burning wave. Some variants of possible stability 
loss due to so-called blow-up modes, or in other words, anomalous nuclear fuel temperature 
and neutron flow evolution are discussed and are found to possibly become a reason for 
trivial violation of internal safety of traveling wave reactor. 

1 Introduction 

Despite the obvious and unique effectiveness of nuclear energy of the new generation there 
are difficulties of its understanding caused by the nontrivial properties of an ideal nuclear 
reactor of the future. 

First, nuclear fuel should be natural, i.e. non-enriched uranium or thorium. Second, 
traditional control rods should be absolutely absent in reactor active zone control system. 
Third, despite the absence of control rods, the reactor must have the property of so-called 
internal safety. This means that under any circumstances the reactor active zone must stay 
at a critical state, i.e. sustain a normal operation mode automatically, with no operator 
actions, through physical causes and laws, that naturally prevent the explosion-type chain 
reaction. Figuratively speaking, the reactors with internal safety are "the nuclear devices 
that never explode" [1]. 

Surprisingly reactors that meet such unusual requirements are really possible. First 
the idea of such self-regulating fast reactor in general form (the so-called breed-and-burn 
mode) was expressed by Russian Feynberg and Kunegin during II Geneva conference 
of 1958 [2] and relatively not long ago "reanimated" in a form of self-regulating fast 
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reactor in traveling nuclear burning wave mode by Russian physicist Feoktistov [3] and 
independently by American physicists Teller, Ishikawa and Wood [4]. 

The main idea of the reactor with internal safety is that fuel components are chosen 
in a way, first, to make the characteristic time Tp of nuclear burning process of active fuel 
component (the fissile component) significantly larger than the time of delayed neutrons 
appearance. And second, in its operation mode all self-regulation conditions should be 
sustained. This is always true if equilibrium concentration hfis of active fuel compo- 
nent, according to Feoktistov wave mode existence condition, is greater than its critical 
concentration^ ?^c^^^ [3]- These conditions are very important, though they are almost al- 
ways practically implementable in case the chain of nuclear transformations of Feoktistov 
uranium-plutonium cycle type [3] is substantial among other reactions in the reactor: 

or the Teller-Ishikawa-Wood thorium-uranium cycle type [4] 

232T/l(n, 7) ^ 233^;^ ^ 233p^ ^ 233 fission), (2) 

In these cases the fissionable isotopes form (^^'^Pu in (1) or ^33^ ^2)) that are the 
active components of nuclear fuel. Characteristic time of such reaction is the time of the 
corresponding /3-decays, that approximately equals to rg = 2.3/ln2 ?a 3.3 days in case (1) 
and ~ 39.5 days in case (2) which is many orders higher than the corresponding time 
for delayed neutrons. 

Self-regulation effect of nuclear burning process is provided by the fact that the system, 
being left by itself, cannot surpass the critical state and enter an uncontrolled reactor 
runaway mode, because the critical concentration is limited from above by finite value of 
equilibrium concentration of active fuel component (plutonium in (1) or uranium in (2)): 
hfis > ricrit (condition of Feoktistov wave mode existence [3]). 

Phenomenologically the process of nuclear burning self-regulation effect is as follows. 
Any increase in neutron fiow leads to quick active fuel component (plutonium in (1) or 
uranium in (2)) nuclei burn-out, i.e. to reduction of their concentration and consequently 
to reduction of neutron flow, though formation of new nuclei by the corresponding active 
fuel component proceeds with prior rate during rg time. On the other hand, if neutron 
flow drops due to external impact, the burn-out speed reduces and consequently active 
fuel component nuclei generation rate will increase followed by increase of a number of 
neutrons generated in the reactor during approximately the same time Tp. 

The details of system of kinetic equations for nuclei (components of nuclear fuel) and 
neutrons (in diffuse approximation) in such chains are rather simple and differ only in the 
depth of description [3-29] of all the possible spectrum of active components of nuclear 
fuel and non-burnable poison^. Fig. 1 with no loss of generality presents characteris- 
tic solutions for such problem (equations (3)- (9) in [29]) in a form of soliton-like waves 
of nuclear fuel components and neutrons concentration for uranium-plutonium cycle in 
cylindrical geometry case. Within the theory of soliton-like fast reactors it is easy to show 
in general case that phase velocity u of soliton-like neutron wave of nuclear burning is 
defined by the following approximate equality [29]: 

""'^-^ ' ^'expf-^a^Y a2 = ^- , , (3) 



2L V^vvr/ V 9^ / 4 hfis-ricrit 



^Concentrations of active element {^^^Pu and ^33^ cycles (1) and (2)), are called equilibrium or critical 
when during the nuclear cycle an equal number of active element nuclei or neutrons is born and destroyed at 
the same time respectively. 

^Here by poison we mean the oxygen nuclei or other elements, chemically bound to heavy nuclides, con- 
struction materials, coolant and the poison itself, i.e. nuclei added to initial reactor composition to control 
neutron balance. 
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Figure 1: Concentration kinetics of neutrons, ^^^U, ^'^^U, ^^^Pu in the core of cylindrical reactor 
with radius of 125 cm and 1000 cm long at the time of 240 days. Here r is transverse spatial 
coordinate axis (cylinder radius), z is longitudinal spatial coordinate axis (cylinder length). 
Temporal step of numerical calculations is 0.1 s. Adopted from [29] 



where hfis and Ucrit are the equilibrium and critical concentration of active (fissile) isotope, 
L is the mean neutron diffusion length, - delay time, associated with the birth of 
active (fissile) isotope and equal to an effective /3-decay period of intermediate nuclei in 
uranium-plutonium cycle of Feoktistov type (1) or in thorium- uranium cycle of Teller- 
Ishikawa-Wood type (2). 

Let us note that expression (3) automatically incorporates a condition of nuclear burn- 
ing process self-regulation, since the fact of a wave existence is obviously predetermined 
by inequality fifis > Ucrit- In other words, the expression (3) is a necessary physical 
condition of soliton-like neutron wave of nuclear burning existence. Let us note for com- 
parison that the maximal value of phase velocity of nuclear burning wave as follows from 
(3) is characterized both for uranium and thorium by the equal average diffusion length 
(L ~ 5cm) of fast neutrons (1 MeV) and is 3.70 cm /day for uranium-plutonium cycle (4) 
and 0.31 cm /day for thorium-uranium cycle (2). 

Generalizing the results of a wide number of numerical experiments [5-9, 12, 14, 16, 18- 
23, 26-29], we can positively affirm that within theory of self-regulating fast reactor in 
traveling wave of nuclear burning mode or in other words, the traveling wave reactor 
(TWR), the principal possibility of main stationary wave parameters control, i.e. its speed 
or thermal power and final fluence both to increase and to decrease them, is determined. 
Obviously it is achieved, according to (3), by varying equilibrium and critical concentration 
of active fuel component, i.e. by purposeful change of the initial nuclear fuel composition. 

In the meantime, the technology problems of TWR are actively discussed in science. 
The essence of these problems on one hand comes to a principal impossibility of such 
project realization, and on the other hand is defined by the following insuperable flaws: 

• High degree of nuclear fuel burn-up (over 20% in average) leading to the following 
adverse consequences: 

— High damaging dose of fast neutrons impact at constructional material (~500 dpa)^; 

^For comparison ~ in Russia for FN-800 the project parameters are stated as 93 dpa. On the other hand it 
is known that one of the main tasks of radioactive stabihty of construction materials investigations in Bochvar 
Hi-tech Institute for non-organic materials (Moscow), is dedicated to achieve 133 to 164 dpa before 2020! 
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— High gas release, requiring increase of gas void length inside excessively long 
fuel rod. 

• Long active zone defining the use of long fuel rods, that have unacceptable param- 
eters from technological use point of view. For instance, this relates to parameters 
characterized by significant increase in: 

— Value of a void effect of reactivity; 

— Hydraulic resistance; 

— Energy consumption for circulation of coolant through the reactor. 

• The problem of nuclear waste due to unburned plutonium processing and nuclear 
waste utilization. 

The main goal of this paper is to solve the specified technological problems of TWR 
on the basis of technical concept enabling, in particular, to increase the damaging dose 
of fast neutrons for in-reactor devices (first of all - fuel rods jackets, reflection shield 
and reactor pit) over ~ 200 dpa. Herein the main technical concept is to achieve the 
given neutron flux on in-reactor devices by defining the speed of fuel movement relative 
to nuclear burning wave speed, the values of which (neutron flux, wave speed and fuel 
movement speed) are in their turn predetermined by choosing the parameters (equilibrium 
and critical concentrations of active component) of initial nuclear fuel composition. 

The article content is following. Chapter 1 is dedicated to brief analysis of theoretical 
development level of self-regulating fast reactor in traveling burning wave mode idea. On 
its basis the problem is formulated and ways to solve it are marked. Chapter 2 considers 
analytical solution for non-stationary ID reacor equation in one-group approximation 
with negative reactivity feedback or, in other words, ID model of Van Dam [7] , yielding 
expression not only for amplitude (pm, phase a and phase speed u, but also for dispersion 
(FWHM) of soliton-like burning wave. Knowing FWHM we may further estimate spatial 
distribution of neutron flux and therefore a final neutron fluence. Chapter 3 is dedicated to 
description of nontrivial dependence of neutron fluence on phase velocity of solitary burn- 
up waves in case of fissible and non-fissible absorbents revealing a possibility of purposeful 
(from the point of view of required neutron fluence and nuclear burning wave speed 
values) change in initial nuclear fuel composition. Chapter 4 analyses the dependence 
of damaging dose on neutron fluence, phase velocity and dispersion of solitary burn-up 
waves. Chapter 5 considers possible causes of internal safety TWR violation, caused by 
Fukushima plutonium's effect, or, in other words, the temperature blow-up modes driven 
either by temperature or neutron flux. Chapter 6 is dedicated to analysis of practical 
examples of the temperature blow-up modes in neutron-multiplying media, and also idea 
of impulse thermo-nuclear TWR is proposed. The conclusion of the paper is presented in 
Chapter 7. 

2 On entropy and dispersion of solitary burn-up 
waves. 

Therefore, below under the specified point of view we discuss physical causes, defining 
main characteristics of soliton-like propagation of " criticality" wave in initially undercrit- 
ical environment, characterized by infinite multiplication factor koo less than unit. It is 
obvious that supercritical area {koo > 1) must be created by some sort of external neutron 
source, e.g. by an accelerator or another super-critical area. In general case, supercritical 
area is a result of breeding effects in fast nuclear systems or burning of fissible absorbents 
(fuel components) in thermal nuclear systems^. Due to progressive burn-out of neutron- 
multiplying environment in supercritical area the area looses its supercritical properties 

^W. Seifritz (1995) first found theoretically a nuclear solitary burn-up wave in opaque neutron absorbers [30]. 
The works of Akhiezer A.I. and his collaborators are dedicated to detection and analysis of super-criticality 
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because koo becomes less than unit. Usually in this case the wave should have stopped 
and diminish. But because neutrons, appearing in source area during breeding diffusely 
"infect" nearby areas, the "virgin" area before the wavefront is forced to obtain the prop- 
erties of super-criticality and the wave moves forward this direction. It is obvious that 
for stable movement of such soliton-like wave there should be some kind of stabilizing 
mechanism, e.g. negative self-catalysis or any other mechanism of negative feedback^. 
In traditional nuclear reactors it is so-called mechanism of negative reactivity feedback. 
Therefore we shall discuss such an example below. 

Let us writhe a non-stationary ID equation of reactor in one-group approximation 
[35-37] with negative reactivity feedback: 

where 93 is the neutron flux [cm~^.s^^]; D - diffusion coefficient, [cm]; jip - reactivity, 
dimensionless value; - total macroscopic absorption cross-section, [cm~^]; v - neutron 
speed, [cm • s~^\. In this case the negative feedback 7 is defined mainly by the fact 
that (4) uses the neutron infinite multiplication factor and therefore flux density must be 
corrected. Let's search for the solution in autowave form: 



ip{x,i) = ip[x-ut) = Lp{i), (5) 

where u is the wave phase velocity, ^ is the coordinate in a coordinate system, which 
moves with phase speed. In such case: 

V dt V dt ^ ' 

As is known [7], the relation u/v by order of magnitude equals to 10^^^ and 10^^^ 
for fast and thermal nuclear systems respectively. Therefore it is obvious that partial 
derivative with respect to time in (4) may be neglected without loss of generality. Further 
taking into account (5) the equation (4) may be presented in this form: 



i'T^ + [fcooW- 1+7^] 9^ = 0, (7) 



where L = {D/Tia)^^"^ is the neutron diffusion length, and is a so-called neutron fluence 
function: 



'4^{x,t) = J (p{x,t')dt'. (8) 


In order to find some physically comprehensible analytic solution of (7) by substituting 
(8) we need to define some realistic view of function A;oo('0)) that is called a burn-up 
function in literature. As is known, the real burn- up function A:oo('0) has a form of some 
asymmetric bell-shaped dependency from fluence value, normalized to its maximal value 
tpmax (fig- 1), therefore, following [7], let us define it without lose of generality in a form 
of parabolic dependency: 

koo — ^max ~l~ (^0 kmax) i ~i 1 | i (9) 

\ Wmax J 



wave in thermal nuclear reactors [31-34], that were first to show the possibility of existence both fast [31-33] 
and slow [34] modes for nuclear burning distribution, or, in other words, super-criticality waves using diffuse 
approximation. 

According to Van Dam [7], the procedure of introduction of reactivity to the non-stationary equation of 
the reactor in one-group approximation though in impUcit form but approximately considers kinetic equation 
of burn-out, particularly plutonium generation in U-Pu cycle or uranium in Th-U cycle. 
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Figure 2: Asymmetric burn-up function as characteristic for realistic burn-up function. Adapted 
from [7]. 



where kmax and kg - are maximal and initial neutron multiplication factors. Substituting 
(9) into (7) we obtain: 



C 



(10) 



where Pmax = kmax - 1, (5 = kmax - ko, 7o = 7- 

Suppose it is necessary to find a partial solution of (10) and write it in the following 
form: 



Pmax + 7oy^ -S 



ip = 0. 



Let us introduce a new unknown function: 

+0O 



/ J 



{ID 



that due to its non-negativity at interval ^ G [0,oo], must satisfy the following boundary 
conditions: = for = 0, oo. 
The equation will take a form: 



— r Pn 



de 



70 — I — 1 



dC 



di 



0. 



(13) 



In order to find a partial solution of (13), we require that the following additional 
condition must be true: 
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Pma.-10^^-S(^-1] =/(x(0), (14) 



where /(^) - yet arbitrary function, the form of which will be defined later. The 
condition (14) is chosen because if it is true, then equation (13) may be integrated. Really, 
if (14) is true, the equation (5) takes the following form: 



i^'^ + ZWOlf^-o, (15) 

That allows us reduce the order of the equation: 



de 



L'^^ + F^{xm = C. (16) 



Here Fi{x) denotes a primitive of f{x) ■, and C is an arbitrary integration constant. 
The order of (16) may be further reduced by multiplying both sides of the equation by 
dx{i)/di: 

f (^) ' + " = ^' ^^^^ 

F2{x) here denotes a primitive of -Fi(x), i-e. "the second primitive" of the function 
f{x) introduced into (5) , and B - a new integration constant. 

The full equation (17) is a separable equation and may be rewritten in the following 
form: 

= ±^= '^^ (18) 



^{B-F2 ixiO) + CxiO) 

On the other hand, (14) may also be considered as a separable equation relative to 
x(^). Then, separating variables in (14) we obtain: 

d^ = . (19) 

■yn (. mb^ 1^ -7(1 



70 70 \Utpm. J 70 

As far as (18) and (19) are equations for the same function x{Cj- Comparing them we 
conclude that the following relation must be true: 



± \l^,{B-FMi)) + Cx{i)) = ^ _ A _ iV _ 1/ (^(^)) . (20) 

V -L^ 70 70 \uipm J 70 

In order to simplify (18) and (19) let's define /(x) in polynomial form of x- The order 
of this polynomial is n. Then i^2(C)> obtained by double integration of f{x) will be a 
polynomial of order (n + 2). Taking a square root, according to (20) should also lead to 
polynomial of order n. therefore n + 2 = 2n=>n = 2. 

Consequently, with the prior presumptions the function /(x) may only be a second- 
order polynomial: 

/(X) = S2X^ + siX + so, (21) 

where sq, si, S2 are polynomial coefficients. 
Double integration of (21) leads to: 

F2{X) = ^X' + 'ix' + ?X' + cix + c2, (22) 
12 2 
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where ci and C2 are the integration constants. Substituting (21) and (22) into (20) we 
get: 



1) - — (s2X^ + siX + so) 



70 70 yu-ilJm J 70 
2 



L2 \B - [-^X" + -^X' + -fx' + ciX + C2)+Cx) (23) 
Further in (23) we set the coefficients at the same orders of x equal: 

+ ^^V = -^. (24) 



-fou'^lpm 70 / ' 

' +'A(r^-'A=-4k. (25) 



lOU^i^m 10 J XlOUi^m 70/ 3L2 ' 

+ — — -— =-T9; (26) 



7o'"V'm 70/ VTo^^^m 70/ V 70 70 70/ L'^ 

2f^--)('— ----)=-('w + w). (27) 
V7owV'm 70 / V 70 70 70 / V ^ / 

^-^-^V = ^S-^. (28) 
70 70 70 / L'^ 

Let's note that first three equations are enough to find the coefficients sQ) sij •52, and 
the remaining two equations may be satisfied by appropriate constants B, C, ci, C2- 

So = Pmax - S; (29) 

25 70 



si 



(30) 



«2 - 777717 /r „ - 777 - 727;^- l-^lJ 



jho 7g^ 

SLnV'mV'^ - Pmax u^V'^ 

After finding sq, si, S2 from this equation system we may consider equation (19) in 
detail, that with accepted notation takes the form: 

= X 7 — X 7 (32) 

& _(_ £2 ] -^2 _|_ ( 23 I £1 ) Y I ( pmax 5 £0 ' 

^ 70Ji^Vm 70 / V 70"V'm 70 / ^ V 70 70 70 , 



Solving this equation yields: 



(X - KY - M2 

where 



-Ndi, (33) 



^ _ ■yputpm 270 /„^N 
A _|_ S2 ' ^ ^ 



Tcru^V^ ' 70 

pmax 5 SO 



\ 70«2V'2, 70 / I 70^2^2^ + 70 
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N= ^ + -. (36) 

Let us introduce a new variable xi into (33) by substituting: 

X-K = Mxi, dx = Mdxi. (37) 
Then, equation (33) will take the following form: 

dxi 
{Xi? - 1 

Hence 



-MNdi. (38) 



Xi = - tanh {MNS, - D) , (39) 
where D is the integration constant. Taking into account x~ K = Mxi- 

X = K- M tanh {MN^ - D) . (40) 
Considering (12), we obtain soliton-like solution in the form: 

ip{C} = M'^Nsech'^iMN^- D). (41) 

Let us remind that when we introduced earlier a new unknown function x{C) (see (12)) 
we obtained an obvious condition for this function as well: 

hm x(0 = 0. (42) 

Let us show, that finally this condition leads to condition of autowave existence. It is 
obvious, that condition (42) considering (40) 



leads to 



lim x(0 = lim [K - Mtanh {MN^ - Z?)] = (43) 



K = M. (44) 



This relation enables us to define values for amplitude 99m, phase a and phase velocity 
u of soliton-like wave: 

a = MN= -HEI^. (45) 

ipm = M^N = ^ ~ ^^'""^ = ~ ^ - (46) 

270 2|7o| 

u = ——. 47 

Prom the condition of nonnegative width (45) and amplitude (46) of nuclear burning 
wave, follows the condition of ID autowave existence, or in other words so called ignition 
condition by van Dam [7]: 

3pmax - 5 = 2kmax + /cq - 3 ^ 0, where 1 - ko > 0. (48) 

We should note here that first analogous results for nonlinear one-group diffusion ID- 
model (4) with explicit feedback and burn-up effects case were obtained in Van Dam 
paper [7]. The same results (see (45)-(47)) were obtained by Chen and Maschek [12] 
while investigating 3D-model by Van Dam using perturbation method, with the difference 
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that the value of neutron fluence associated with the maximum of burn-up parameter 
koo adapted to the transverse buckhng. In other words they considered the transverse 
geometric buckhng mode as a basis for perturbation. Hence they introduced a geometric 
multiphcation factor kcB due to transverse buckhng into two-dimensional equation (4), 
that led to change in some initial parameters {pmax — kmax — kQB^ S — kmax — 

ko) and 

therefore led to change of conditions of autowave existence in 3D case: 

Spmax - S = 2kmax + ko - SkcB > 0, where kcB - A;o > 0, (49) 

that in case kQB = 1 is exactly the same as the so-called ignition condition by Van Dam [7] . 

Obviously from the point of view of the more detailed Feoktistov model [3] analysis, 
thoroughly considered in [29] and connected to concept of internal safety of nuclear sys- 
tems, the condition (49), on the one hand is a necessary but non sufficient, and on the 
other hand is a hidden form of necessary condition of wave existence according to Feoktis- 
tov, where equilibrium concentration fifia of active fuel component must be greater than 
its critical concentration ricrit {"^fis > ncrit) [3, 29]. Physics of such hidden at first sight 
but simple interrelation will be discussed later (see Chapter 3). 

Returning to ID reactor equation solution (10) in one-group approximation with neg- 
ative reactivity feedback, let us write it in more convenient form for analysis 

'P = fm sec /i^(a^) = (pm sec [a{x — ut)]. (50) 

where ifm - amplitude of neutron flux; 1/a - characteristic length proportional to the 
soliton wave width, that is a full width at half-maximum (FWHM) by definition and 
equals to 



Ay2 = FWHM = 2ln(l + V2ja-^, [cm]. (51) 
It is obvious that integrating (49) yields area under such soliton: 

Aarea = 2ipm/a, [cm^]. (52) 

In order to estimate the extent of the found parameters infiuence on the dynamics 
of a soliton-like nuclear burning wave stability let us invoke an information-probability 
approach, developed in paper by Seifritz [38]. Hence, following [38], let us write down 
the expression for mean value of information or, more precisely, entropy of the considered 
process: 



S = —kB J p{x) In p{x)dx, (53) 

— oo 

where p{x) is the function of probability density relative to dimensionless variable x; 
lnl/p{x) - mean information value; kB - Boltzmann constant. 
Substituting soliton-like solution (49) into (53) we obtain 

oo oo 

S = -2kB I sec h\a£,) In [sec h\aO] d ( —) = AkB ! ^"^^"^^^dy, (54) 
J \ur^J J cosh"=(/3y) 



where (3 = aluTp) is another (dimensionless) scaling factor, t/j - proper /3-decay time 
of the active component of nuclear fuel. Note the procedure of making the x argument 
dimensionless (54), which considers the fact that the amplitude of neutron fiux is propor- 
tional to the phase velocity of nuclear burning considering (47), i.e. ~ u. Calculating 
integral (54) leads to the following quite simple expression for the entropy 

, 4(1-1..2) , fa 4Ml-ln2)2L 

/3 auTp y/kcB - ko UTp 
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Figure 3: (a) Weakly dispersive wave pattern, obtained by Chen and Maschek [12] investigating 
a 3D-model by van Dam using perturbation method. Example with the following parameters: 
L = 0.02 m, kcB = 1-04, ko = 1.02 and k^ax = 1-06; 0„ = 10^^ m~^s~^, u = 0.244 cm/day; (b) 
Strongly dispersive wave pattern, obtained Chen et. al [18] within a 3D-model of traveling wave 
reactor. Example with the following parameters: L = 0.017 m, kcB = 1.00030, fco = 0.99955; 
= 3 • 10^^ m~^s~^, u = 0.05 cm/day. 



that in case 

S ~ kB = const, (56) 

UT/s 

points to isoentropic transport of nuclear burning wave. 

It is interesting to note here that if the width A1/2 0, then due to isoentropicity of 
the process (56), the form of soliton becomes similar to so called Dirac J- function. 

If we introduce two characteristic sizes or two length scales {h = 1/a and I2 = (pm), 
then we may see, that when the first of them is small, the second increases and vice 
versa. This happens, because the area Aarea (45) under soliton must remain constant 
because Aarea oc lil2- In this case the soliton entropy approaches zero, because entropy 
is proportional to relation of these values {S cx h/h — )• 0). Denoted properties are the 
consequence of the fact, that the scale h is a characteristic of dispersion of the process, 
while another scale I2 is a characteristic of soliton non-linearity. If li <^ I2, then the 
process is weak-dispersion (fig. 3a), and vice versa, if h ^ h, then the process is strong- 
dispersion (fig. 3b). In the last case the soliton amplitude becomes relatively large (case 
of 5-function). And finally, if h = I2, then soliton speed u oc (^i^i)"*^^^ (^^^ (48)) is 
proportional to geometrical mean of dispersion and non-linearity parameters. 

On the other hand, it is obvious, that according to (45) and (47) the burning wave 
width is a parameter that participates in formation of the time for neutron fluence 
gathering Tp on the internal surface of cladding material of TWR long fuel rod 

~ 'Pm = ■ ^m- (57) 

U 

And finally, one more important conclusion. From the analysis of (45) and (46) it 
is clear that initial parameter ko for burning zone (fig. 4) is predefined solely by nuclear 
burning wave burn-up conditions, i.e. by parameters of external neutron source and burn- 
up area composition (fig. 4). In other words, this means that by tuning the corresponding 
burn-up conditions for a given nuclear fuel composition we can set the certain value 
of nuclear burning front width. Moreover, by selecting the corresponding equilibrium 
hfis and critical ncrit concentrations of active nuclear fuel component, we can define the 
required value of nuclear burning wave speed u. It becomes obvious how we can control 
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Figure 4: Schematic sketch of the two-zone cyhndrical TWR. 



the corresponding neutron fluence gathering time on the internal surface of cladding 
material of a TWR long fuel rod. 

Therefore we may make an important conclusion that for realization of a real TWR 
project that would really have the property of internal safety we need knowledge about the 
physics of nuclear burning wave burn-up and a knowledge about the interrelation between 
the speed of nuclear burning wave and the fuel composition, the properties of which, as 
shown in [29], are completely defined by the values of the equilibrium hfis and critical 
ncrit (see (3)) concentrations of the active nuclear fuel component. Below we consider this 
in detail. 



3 Control parameter and condition of existence 
of stationary wave of nuclear burning. 



According to the information given above a question arises: "What does nuclear burning 
wave speed in uranium- plutonium (1) and thorium- uranium (2) cycles mainly depend 
on?" The answer is simple and obvious. The nuclear burning wave speed in both cycles 
(far away from burn- up source) is completely characterized by values of equilibrium njjs 
and critical ricrit concentrations of active nuclear fuel component. 

First of all this is determined by a significant factor that value of equilibrium h fig and 
critical ricrit concentrations of active nuclear fuel component completely identify neutron- 
multiplying properties of the fuel environment because they are the conjugate pair of 
integral parameters, that due to physical content fully and adequately characterize all 
the physics of nuclear transformations, predefined by initial fuel composition. This is also 
simple to see from a simple analysis of kinetics equations solutions for neutrons and nuclei, 
used in different models [5-9, 12, 14, 16, 18-23, 26-29]. Hence the conclusion is following. 
Independently of nuclear cycle type and initial fuel composition the nuclear burning wave 
speed is defined by equilibrium hfis and critical ricrit concentrations of active nuclear 
fuel component through the so-called para-parameter a (see (3)). And, consequently, as 
numerical experiments results show, (fig. 5 [29]), follows Wigner statistics. 

This is corroborated by the following fac. In paper [20, 21] a boundary-value problem 
was investigated for stationarity of nuclear burning wave, formulated within diffusion 
equation for fluence of neutrons and equations for kinetics of nuclear density of mother 
and daughter nuclides: 






(59) 
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where 

l/'end 

G(iV,v) = I g[N)dil^ (60) 



is a function^ of neutron fluence generation, g{N) - function of neutron generation, that is 
a Unear function of nuclei concentration Ni, N - column of mother and daughter nuclides 
Ni, ij) - neutron fluence, ipend ~ some final neutron fluence, v - average neutron speed, 
D - diffusion coefficient, a -matrix of microscopic neutron absorption and capture cross- 
sections, and A - matrix of radioactive decay constants of /3-active nuclei. 

It is shown [20], that boundary- value problem for stationary nuclear burning wave 
for equations (58)- (59) is a spectral non- linear differential problem. As a spectral (free) 
parameter, a dimensionless speed W of nuclear burning wave was selected. Such presen- 
tation of a problem allows analytic investigation of the final fluence change behavior for 
a nuclear burning wave due to change of absorbent properties, concentration of which 
according to [20, 21, 26, 27], controls (in a zero approximation of perturbation theory) 
the value of dimensionless wave speed: 

^^UT2^1 atW^l, (61) 

L b 

where u -nuclear burning wave speed, T2 - internal time scale, equal to characteristic 
time of /3-decay of intermediate nuclide {t2 = 3.47 days for "^^^Np in U-Pu cycle and 
T2 = 36.6 days for ^^^Pa in Th-U cycle), L is the neutron diffusion length, 1/6 is the 
linear dependency slope (61), j? - dimensionless effective concentration of absorbent, pQ - 
upper limit of absorbent concentration area where the wave can exist. 

Here we obtain an important result showing that the final fluence and absorbent 
concentration in extremal case W ^ (zero order of perturbation theory) are the solutions 
of two equilibrium conditions for a stationary nuclear burning wave: 

^The expression (59) should be understood as integral along system path in configuration space of variables 
{N,ip) for the given spatial point 



13 



gd^ = 0, (62) 



Ipend 

M{lPend)= j (Vend - V) ffdV = 0. (63) 


For more clear understanding of the physical meaning of these conditions, the authors 
of [20] propose a simple, but elegant and deep analogy, that the obtained conditions exactly 
coincide by form and accuracy with the conditions for a lever subject to a perpendicular 
distributed force ^(V)! that is applied on a segment ^ V ^ Vend along the lever. In other 
words, in this analogy, the conditions (62) and (63) represent conditions of zero force sum 
and momentum sum respectively. Therefore, if the first condition is an integral condition 
of the neutrons generation and absorption equality , then, according to the analogy, the 
second condition may be called the condition of the neutrons generation and absorption 
"momenta" equality. 

If expression for g ^iV^ is presented in the form of a sum of the contributions of fuel 
g^ and absorbent g"^ = p in neutron generation: 



9 = 9^-9^, (64) 

iccording to conditions (62) 
le mean "momenta" of neut 

form [26] 



then according to conditions (62) and (63) it is easy to obtain a modified condition 
(63) for the mean "momenta" of neutrons generation (^M^'j and absorption (^M^^ in the 



f g^Tpdij) f g^il^dij) 

(M^ (Vend)) = = ^ = (M^ (Vend)) , (65) 

Wend Wend 

J g^d^P J g^di; 



which by definition (62) has a trivial root V = Vend- 

The main advantage of such presentation of equation (65) is the notable fact that its 
left part, i.e. (M^), is defined solely by properties of the fuel, while its right part, i.e. 
^M"^), is defined solely by absorbent. Moreover, according to [26], the initial absorbent 
concentration is absent in the equation (65). The plots for left and right parts of equation 
(65) for different cases are presented at fig. 6 [26]. 

Simple analysis of fig. 6 shows that with increase of concentration of non-burnable 
absorbent (curves 2-5 at fig. 6a) the value of the final Vend fluence increases, but, according 
to (61), nuclear burning wave speed decreases. And vice versa, with increase of non- 
burnable absorbent concentration (curves 1-4 at fig. 9b) the value of the final fiuence Vend 
decreases, while the velocity of nuclear burning wave at p = also decreases, but in a 
complex manner (see fig. 7), depending on the slope behavior (61). Increase in absorbent 
concentration (curves 5-6 at fig. 6b) does not change the effect qualitatively, but leads 
to a strong quantitative change. Here it turnes out that the final fluence is very sensible 
to concentration of burnable absorbent: curve 6 at fig. 6b for burnable absorbent crosses 
curves 1-4 at much greater values of final fluence. 

Therefore, in paper [26, 27] on the basis of results of numerical solution of boundary- 
value problem for stationary nuclear burning wave and developed earlier analytical theory 
of nuclear burning wave [20, 21] it is shown (in a zero order of perturbation theory ac- 
cording to W), that control parameter, enabling us to change speed of nuclear burning, 
and hence the reactor power, is the effective concentration of absorbent. Moreover, this 
control parameter may both increase and decrease the specifled parameters. According 
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Figure 6: Graphical solution of equation (65) for final fluence ipend with (a) different speeds 
of absorbent burn-out and (b) different enrichment of "^^^Pu for non-burnable and burnable 
adsorbents. At fig. 7a curve 1 is a graph of right equation part, and curves 2-5 - left equation 
part for the absorbent cross-sections 0, 1, 2 and 3 barn respectively. At fig. 7b curves 1-4 are 
plots of right part for initial concentrations of ^^^Pu 0, 3, 5, 7% respectively, and curves 5 and 
6 - left equation part for the absorbent cross-sections and 2 barn respectively. Adopted from 
[26]. 
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Figure 7: Dependency of nuclear burning wave speed on concentration of non-burnable ab- 
sorbent p, (a) considering burning-out of ^^^Np (curve 1) and not considering burning-out of 
A^p (curve 2), (b) for different initial fuel compositions: 1 - 100% ^^^U; 2 - 93% ^^^U + 7% 
235^7; 3 - 93% [/ + 7% 239p^^ Adopted from [26]. 
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to [20, 21, 26, 27] ". . .this is achieved by purposeful change of initial composition of the 
reactor" . 

This conclusion of authors [20, 21, 26, 27] requires the following comment. Due to our 
earnest conviction the results of these works are doubtlessly very interesting and informa- 
tive, but as strange as it may seem, the conclusions may not be considered satisfactory. 
And this is why. 

Really, as shown in [29] , the control parameter is not effective absorbent concentration 
(p) (or nuclear burning wave speed (W), or maximal neutron fluence (V'end))) but a so- 
called para-parameter of nuclear TWR-burning that forms as shown above (see (3)), by a 
conjugated pair of integral parameters or, in other words, by the equilibrium and critical 
concentrations of active nuclear fuel component. It is important to note here that despite 
the fact that each of this parameters pair are variable during nuclear burning, but their 
ratio 



is not only a characteristic constant value for the given nuclear burning process [29], that 
due to its physical meaning fully and adequately characterizes all the physics of nuclear 
transformations predefined by initial fuel composition, but also determines, first (and it is 
very important!) the conditions of nuclear burning wave existence (3), second, the neutron 
nuclear burning wave speed (see (3)) and, third dimensionless width (66) of super-critical 
area in the burning wave of nuclear fuel active component. 

Therefore, when authors [20, 21, 26, 27] state that they control the values of the 
parameters with purposeful variation of initial reactor composition, it really means that 
changing the effective concentration of absorbent, they purposefully and unambiguously 
vary (by definition (see [3, 29])) the equilibrium and critical concentrations of nuclear 
fuel active component, or the control para-parameter (66) of the nuclear TWR-burning 
process ( which means the same). 

Here we should note that the adequate understanding of the problem of control pa- 
rameter discrimination and determination is not a simple or even scholastic task, but ex- 
tremely important for effective solution of another no- less-important problem (in essence, 
the main one!), related to investigation of the nuclear burning wave stability conditions 
problem. In our opinion, the specified conditions for stationary nuclear burning wave 
existence (3) and (62)-(63) obtained respectively in papers [29] and [20, 21, 26], reveal 
the path to a purposeful application of the so-called direct Lyapunov method [39] being 
the base theory for movement stability, and therefore, the path to a reliable justification 
of the Liapunov functional minimum existence (if it exists!) [39-42]. Some variants of a 
probable loss of solution stability due to anomalous evolution of nuclear fuel temperature 
are considered in Chapter 5. 

Meanwhile, despite all said above we can draw a conclusion that "differential" [29] 
and "integral" [20, 21, 26] conditions for stationary nuclear burning wave existence are 
not only a reliable source for complete description of the wave reactor physics , but in the 
future can become a basis for engineering calculations of contemporary project of TWR 
with optimal or preassigned wave properties. 



4 On dependence of damaging dose from neutron 
fluence, phase velocity and dispersion of solitary 
burn- up waves 



As follows from expression for soliton-like solution (50), it is obvious, that it is defined 
by three parameters - maximal neutron flux ^pm, phase a and speed u of nuclear burning 




(66) 
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Figure 8: Swelling of austenitic Phenix cladding compare to ferritic-martensitic materials, ODS 
included. Adopted from [48]. 



wave. Here, even if we assume that we can control them, an obvious question emerges - 
how or by which conditions are the optimal values of these parameters determined? Let's 
try to answer this question in short. 

It is known, that high cost effectiveness and competitiveness of fast reactors, including 
TWR, may be achieved only in case of high nuclear fuel burnup^. Experience of the fast 
reactors operation shows that the main hindrance in achieving high nuclear fuel burnup 
is the insufficient radiation resistance of the fuel rod shell. Therefore the main task of 
the radiative study of materials is not only to study the physics of the process , but also 
to create the material (or select among the existing materials), which would keep the 
required level of performance characteristics being exposed to the neutron irradiation. 
One of the most significant phenomena leading to a premature fuel rods destruction is 
a vacansion swelling of shell material [44-47]. Moreover, absence of satisfactory swelling 
saturation and even its acceleration with increase of damaging dose leads to a significant 
swelling (volume change up to 30% and more) and subsequently to significant increase 
in active zone elements size. Herewith the consequences of such effect are amplified by 
the fact that a high sensitivity of the swelling to temperature and irradiation damaging 
dose leads to distortion of active zone components form because of temperature and 
irradiation gradients. And finally, one more aggravating consequence of high swelling 
is almost complete embrittlement of construction materials at certain level of swelling^. 
Consequently, in order to estimate the possible swelling level, a damaging dose (measured 
in dpa) initiated by fast neutrons, e.g. in fuel rods shell, must be calculated (fig. 8). 

Usually to evaluate the displacements per atom (DPA) created by the spallation 
residues a so-called modified NRT method is applied [49, 50], considering a known Lind- 
hard correction [51]. Within such modified NRT method, the total number of displace- 
ments produced by the residues created by spallation reactions in the energetic window 
can be calculated as the addition of the displacements produced by each of these residues 
{Z,A), that in its turn leads to the following expression for the displacements per atom: 

^Because the maximal burnup of FN-600 reactor is ^ 10% [43], the burnup factor ^ 20% for TWR reactor 
may be considered as ultra reasonable. 

^It is known that fuel rod shell diameter increase due to swelling is combined with anomalously high corrosion 
damage of the shell by the fuel [43] 
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Figure 9: Operating condition for core structural materials in different power reactors [52]. The 
upper yellow inset represents the data of Pukari M. & Wallenius [47]. 



N N 

ndpa = t-Y,{^dpa) J ^{Ei)dEi = Y,{<ydpa)Vit. (67) 

where 

(a^p,) {Ei) = (a, Z, A)) ■ d (E^ (Z, A)) , (68) 

while a(i{Ei, Z, A) is the displacement cross-section of recoil atom {Z,A), produced at 
incident particle energy Ei, ^(E) is the energy-dependent flux of incident particles during 
time t and d {E^{Z, A)) is the number of displacements created at threshold displacement 
energy Ed of recoil atom {Z,A) or its so-called damage function. 

Actually, all the recoil energy of the residue Er is not going to be useful to produce 
displacements because a part of it is lost inelastic scattering with electrons in the medium. 
An estimation of the damage energy of the residue can be calculated using the Lindhard 
factor ^ [51] 

Edam = Eri. (69) 

The number of displacements created by a residue {Z,A) are calculated using this 
damage energy (69) and the NRT formula: 

d{Ed{Z,A))=v^y (70) 

where {Ed) is the average threshold displacement energy of an atom to its lattice site, 
r/ = 0.8 [49, 50]. 

Consequently the condition of maximal damaging dose for cladding materials of fast 
neutron reactors, taking into account metrology data of IAEA [52] (see fig. 9) and con- 
temporary estimates by Pukari M. & Wallenius (see yellow inset at fig. 8) will have the 
form: 
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Figure 10: Discretized displacement cross-section for stainless steel based on a Lindhard model 
and ENDF/B scattering cross- section. Adopted from [53]. 



ridpa ^ {cFdpa) ■ ^ ^ ^ 200 [dpa]. (71) 

It is obvious that in this case the strategy for required wave parameters and allowable 
values of neutron fluence selection for future TWR project must take into account the 
condition for maximal damaging dose (70) for cladding materials in fast neutron reactors 
and, therefore, must comply with the following dpa-relation: 

{^dpa)-^- — ^ 100 [dpa]. (72) 
u 

Of course, a question arises here whether the parameters of wave and neutron fluence 
that can provide the burnup level of active nuclear fuel component in TWR-type fast 
reactor at least 20% are possible. Since we are further interested in cladding materials 
resistible to fast neutron damaging dose we shall assume that displacement cross-section 
for stainless steel, according to Mascitti et al. [53] for neutrons with average energy 2 MeV 
equals {a dpa) ~ 1000 dpa (fig. 10). 

Analysis for nuclear burning wave parameters in different authors' models of TWR 
[4, 5, 13-15, 18, 25, 26, 28, 29] presented in Table 1, shows that in case of U-Pu cycle 
none of the considered models satisfy the dpa-parameter, while within Th-U cycle two 
models by Teller [4] and Seifritz [5] groups correspond well to major requirements to wave 
reactors. 

On the other hand it is obvious that authors of [13-15, 18, 25, 26, 28] just didn't 
consider or take into account the problem of dpa-parameter in claddings materials, because 
as of the date of these works, they were mainly interested in the fact of wave mode of 
nuclear burning existence in U-Pu and Th-U cycles. 

However, as analysis of Table 1 reveals, the procedure of account for dpa-parameter is 
not problematic, but it leads to unsatisfactory results relating to burn-out of main fissle 
material. 

In other words, from the analysis of Table 1 is follows that when {adpa) ~ 1000 dpa, the 
considered above dpa-condition of maximal possible damaging dose for cladding materials 
of fast neutron reactors 
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Ai/2 

[cm] 


u 

[cm/day] 


[cm~2s~^] 


[cm~2] 


{O'dpa) 

[barn] 


^dpa 

200 


Fuel 
burn-up 


Solution 


U-Pu cycle 


Sekimoto [14] 


90 


0.008 


3.25 ■ 10^^ 


3.2 ■ 10 23 


1000 


3.2 


~43% 


No 


Rusov [29] 


200 


2.77 


10^8 


6.2-1024 


1000 


62 


-60% 


No 


Pavlovich [26] 




0.003 




1.7-1024 


1000 


17 


-30% 


No 


Fomin [15] 


100 


0.07 


2 ■ 10^^ 


2.5-1024 


1000 


25 


-30% 


No 


Fomin [13] 


125 


1.7 


5 ■ 10^' 


3.2-1024 


1000 


32 


-40% 


No 


Chen [18] 


216 


0.046 


3 ■ 10^^ 


1.2-1024 


1000 


12 


-30% 


No 


T. Power [28] 












1.75 


-20% 


No 


Th-U cycle 


Teller [4] 


70 


0.14 


~2 ■ 10^5 


8.6-1022 


1000 


0.96 


-50% 


Yes 


Seifritz [5] 


100 


0.096 


10^5 


9.0-1022 


1000 


0.90 


-30% 


Yes 


Melnik [25] 


100 


0.0055 


0.5-1016 


7.9-1024 


1000 


~80 


-50% 


No 


U-Pu (+ moderator 




Example 


100 


0.234 


2.5 ■ 10^5 


9.2-1023 


100 


0.92 


-20% 


Yes 


Ideal TWR 








1024 


100 


1.0 


^20% 


Yes 



Table 1: Results of numerical experiments on investigations of wave mode parameters based 
on U-Pu and Th-U cycles 



V'looo = — - 10'^ [cm-% (73) 
u 

is not met by any example in Table 1. Here i/'iooo is the neutron fluence in case {o'dpa) ~ 
1000 dpa, if - neutron flux, A1/2 a-J^d u are width and speed of soliton-like nuclear burning 
wave. 

It is obvious that in order to significantly increase burnup level on one hand a neutron 
fluence should be increased by an order, and on the other hand the maximal damaging dose 
for cladding materials should be reduced. Such a controversial condition may be fulfilled 
considered that fuel rod shell construction material radioactive damage level reduction 
for a given amount may be achieved by reducing neutron flux density and energy (see 
fig. 10). The latter is achieved by placing a specially selected substance between fissile 
medium and fuel rod shell, which has some specific characteristics of neutron moderator 
and absorbent. 

On the other hand, within the framework of reactor neutron physics it is known [35, 54] , 
that moderator layer width estimate Rmod is: 

Rmod^^^^--Mp^, (74) 
^5 + ? ^mod 

where « {as) Nmod and {(ja) Nmod neutron macroscopic scattering and ab- 

sorption cross-sections respectively, {as) and {aa) - neutron microscopic scattering and 
absorption cross-sections respectively averaged by energy interval of moderating neu- 
trons from Efuel = 2 MeV to Emod = 0.1 MeV, N^od ~ moderator nuclei density, 
^ = 1 -|- (j4 -|- 1)2 In [{A — 1) /A + 1] /2A - neutron energy decrement of its moderating in 
the moderator- absorbent medium with atomic number A. 

It is clearly seen that the process of neutron moderating from 2.0 MeV to 0.1 MeV 
energy in moderator-absorbent of a given width (see Table 2) creates a new, but satis- 
factory level of maximal possible damaging dose for cladding materials, corresponding to 
{(^dpa) ^ 100 dpa (fig. 10). Therefore if we are satisfied with main fissile material burn-out 
level around —20%, then analyzing Table 1 and Table 2, the conditions accounting for 
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Moderator 


Mass 

number, 

A 


Mean loga- 
rithmic 
energv, £ 


density, p, 
g/cw? 


Impacts 
number 
required for 
moderating, n 


Neutron 
mean free 
path, A 


Moderator 

layer 

width, 

Rraod, cm 


Be 


9 


0.21 


1.85 


11 


1.39 


15.3 


C 


12 


0.158 


1.60 


15 


3.56 


53.4 


H2O 


18 


0.924 


1.0 


2.5 


16.7 


41.6 


H2O + B 








2.5 


10.0 


25.0 


He 


4 


0.425 


0.18 


5.41 


11.2 


60.7 



Table 2: Moderating and absorbing properties of some substances, moderator layer width 
estimate for moderating neutron from Ejuei = 1-0 MeV to Emod = 0.1 MeV. 



dpa-parameter problem and contemporary level of radioactive material science will have 
the following form: 

V'loo = 10'' [cm-% (75) 

u 

where "^loo is the neutron fluence (with 0.1 MeV energy) on cladding materials surface in 
case {(Tdpa) ~ 100 dpa (see fig. 10 and "ideal" case in Table 1). 

And finally we may make the following intermediate conclusion. As shown above in 
Chapter 3, the algorithm for determining the parameters values (73) is mainly determined 
by para-parameter that plays a role of a kind of a response function to all the physics of 
nuclear transformations, predefined by initial fuel composition. It is also very important 
that it unambiguously determines first, the conditions of nuclear burning wave (3), second, 
the neutron nuclear burning wave speed (see (3)) and third, a dimensionless width (66) 
of supercritical area in nuclear fuel active component burning wave. 

Based on para-parameter ideology [29] and Pavlovych group results [20, 21, 26], we 
managed to pick up a mode for nuclear burning wave in U-Pu cycle, having parameters 
shown in the Table 1 satisfying (73). The latter means that the problem of dpa-parameter 
in claddings materials in the TWR-project is currently not an insurmountable technical 
problem and can be successfully solved. 

In our opinion, the main problem of TWR are the temperature blow-up modes that 
take place due to coolant loss as observed during Fukushima nuclear accident. Therefore 
below we shall consider possible causes of TWR internal safety breach due to temperature 
blow-up mode. 

5 Possible causes of TWR internal safety breaches: 
Fukushima plutonium's effect and temperature blow- 
up mode 

It is known that with loss of coolant at three nuclear reactors during Fukushima nuclear 
accident its nuclear fuel melted. It means that temperature in the active zone in some 
moments reached melting point of uranium- oxide fuel^ i.e. ~3000°C. 

Surprisingly in scientific literature today there are absolutely no neither experimen- 
tal nor even theoretically calculated data on behavior of capture cross-sections of ^^^^7 
and ^^^Pu depending on temperature at least in 1000-3000°C range. At the same time 
there are serious reasons to suppose that the cross-section values of the specified elements 

^We should note here, that the third block partially used MOX-fuel, enriched with phitonium 
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increase with temperature increase. At least we may point for example at qualitative es- 
timates of Ukraintsev [55] , Obninsk Institute of Atomic Energetics (Russia) , that confirm 
the possibility for cross-sections growth phenomenon for '^^^Pu in 300-1500°C range. 

Obviously, such anomalous temperature dependency of capture and fission cross- 
sections of ^j^j 239 may drastically change neutron and thermal kinetics of nuclear 
reactors including perspective fast uranium-plutonium new generation reactors (reactors 
of Feoktistov (1) and Teller (2) type), that we classify as fast TWR reactors. Hence it 
is very important not only to know anomalous temperature behavior of '^^^JJ and "^^^Pu 
capture and fission cross-sections, but also it becomes critically important to know their 
influence on heat transfer kinetics, because it may turn into a reason of positive feedback^ 
(PF) with neutron kinetics leading to undesirable solution stability loss (the nuclear burn- 
ing wave) as well as to a trivial reactor runaway with a subsequent nontrivial catastrophe. 

A special case of PF is a non- linear PF, leading to system evolution in so-called blow- 
up mode [56-61] or, in other words, in such a dynamic mode when one or some of modeled 
values (e.g. temperature and neutron flux) grows to infinity at a finite time . In reality 
instead of infinite values a phase transformation is observed in this case, that can in its 
turn become a first stage or a precursor of future technogenic disaster. 

Investigation of temperature dependency of capture and fission cross-sections of 238f/ 
and 239pit in 300-3000° C range and, correspondingly, the kinetics of heat transfer and its 
infiuence on neutron kinetics in TWR is the main goal of the chapter. 

Heat transfer equation for uranium-plutonium fission environment is: 

, ^ dT(f,t) 

p{r,T,t)-c{r,T,t)-—^ = 

= ^ (f, T, t) ■ AT (f, t) + V\i (f, T, t) ■ VT (r, t) + ^ (r, T, t) , (76) 
where the effective substance density is 

p{f,T,t) = Y,Ni{f,T,t)-pi, (77) 

i 

Pi are tabulated values, Ni (r, T, t) are the components concentrations in the medium, 
while the effective specific heat capacity ( accounting for the medium components heat 
capacity values Cj) and fissile material heat conductivity coefficient (accounting for the 
medium components heat conductivity coefficients 'iii{T)) respectively are: 

c (r , T,t) = Y, Ci {T)Ni (f, T, t) , (78) 

i 

i<{f,T,t) = J2MT)Ni{r,T,t). (79) 

i 

Here (r , T, t) is the heat source density generated by nuclear fissions Ni of fissile 
metal composition components that vary in time. 

Theoretical temperature dependency of heat capacity c(T) for metals is known: at low 
temperatures c{t) ~ T^, and at high temperatures c(T) — )• const, notably the constant 
value (const « 6 Cal/{mol ■ deg)) is determined by Dulong-Petit law. At the same time it 
is known that thermal expansion coefficient is small for metals, therefore the specific heat 
capacity at constant volume c^ almost equals to the specific heat capacity at constant 
pressure Cp. On the other hand, the theoretical dependency of heat conductivity Hi(T) at 
high temperature of "fissile" metals is not known, while it is experimentally determined 



^Positive Feedback is a type of feedback when change in output signal leads to the change in input signal, 
that in its turn leads to further deviation of output signal from its original value. In other words PF leads to 
instability and appearance of qualitatively new (often self-oscilation) systems. 
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that heat conductivity coefficient H(T) of fissile medium is a non- linear function of tem- 
perature (e.g. see [62], where heat conductivity coefficient is given for a- uranium 238 and 
for metallic plutonium 239, and also [63]). 

Further while solving thermal conductivity equations we used the following initial and 
boundary conditions: 

T{r, t = 0) = 300 K and jn = H [T{r G t) - Tq] , (80) 

where jn is the normal (to the fissile medium boundary) heat fiux density component, ^(T) 
is the thermal conductivity coefficient, ?R. is the fissile medium boundary, Tq - temperature 
of the medium adjacent to the active zone. 

Obviously, if the cross-sections of some fissile nuclides increase, then due to exothermic 
nuclei fission reaction the direct consequence of significantly non-linear kinetics of parental 
and child nuclides in the nuclear reactor is an autocatalyst increase of generated heat, 
similar to autocatalyst processes for exothermic chemical reactions. In this case thermal 
fiux density qJp (r, <I>, T, t) that characterizes generated heat amount will be: 

4 (f, CD, r, t) = CD (f, T, t) J2 Q{^} (r, T, t) N, (f, T, t) , [W/cnr'] , (81) 

i 

where 

^max 

$ (f, T,t)= j $ (r, E, T, t) dE 



is the full neutron fiux density; cD (f*, E, T, t) - neutron fiux density with energy E; - 
mean generated heat emitted due to fission of one nucleus of the i-th. nuclide; 

a}{r,T,t)= I a}{E,T)p{r,E,T,t)dE 


is the fission cross-section of the i-th nuclide averaged over the neutron spectrum; 

p (f, E, T,t) = ^ (f, E, T, t) (f, T, t) 

is the neutron energy distribution probability density function; cr^fiE, T) is the microscopic 
fission cross-section of the i-th nuclide that, as known, depends on neutron energy and 
fissile medium temperature (Doppler effect [35]); Ni {f,T,t) - density of the i-th. nuclide 
nuclei. 

As follows from (81), in order to build thermal source density function it is necessary 
to solve a problem of theoretical dependency of the cross-sections (f , T, t) averaged over 
the neutron spectrum on the temperature of reactor fuel (fissile medium). As is known, 
the ifnluence of the nuclei thermal motion o in the medium comes to a broadening and 
height reduction of the resonances. By optical analogy this phenomenon is referred to as 
Doppler effect [35]. Since in low energy area the resonance levels are observed only for 
heavy nuclei, then Doppler effect is notable only during the interaction of neutrons with 
such nuclei. And the higher environment temperature the stronger is the effect. 

Therefore a program was developed using Microsoft Fortran Power Station 4.0 (MFPS 
4.0) that allows at the first stage calculate the cross-sections of the resonance neutron 
reactions depending on neutron energy taking into account the Doppler effect. The 
cross-sections dependency on neutron energy for reactor nuclides from ENDF/B-VII 
database [64], corresponding to 300K environment temperature, were taken as the in- 
put data for the calculations. For example, the results for radioactive neutron capture 
cross-sections dependency on neutron energy for '^^^JJ are given in fig. 11 for different tem- 
peratures of the fissile medium in 300K-3000K temperature range. Using this program, 
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Figure 11: Calculated dependency of radioactive neutron capture reaction on its energy for 
ll^U at different temperatures within 300K to 3000K. 



the dependency of scattering, fission and radioactive neutron capture cross-sections for 
essential reactor fuel nuclides q2^U, q2^U, q2^U and q^39Pu for different temperatures in 
range 300K to 3000K were obtained. 

At the second stage a program was developed to obtain the calculated dependency 
of the cross-sections cij (r, T, t) averaged over the neutron spectrum for main reactor 
nuclides and for main reactor nuclear neutron reactions for the specified temperatures. 
The averaging of neutron cross-sections for the Maxwell distribution was performed using 
the following expression: 

{a{Eu^,T)) = , 



where E^^ is the the upper limit of neutrons thermalization, while for the procedure of 
neutron cross-sections averaging over Fermi spectrum the following expression was used: 

oo 

/ a{E,T)E-^dE 

{a{Eii^,T)) = ^""^ ^ , 

/ E'^dE 

During further calculations in our programs we used the results obtained at the first 
stage i.e. the dependency of reaction cross-sections on neutron energy and environment 
temperature (Doppler effect). The neutron spectrum was specified in a combined way - 
below the limit of thermalization En^ the neutron spectrum was described by Maxwell 
spectrum $m (-E'n); above Eum but below Ep (upper limit for Fermi neutron energy 
spectrum) the neutron spectrum was described by Fermi spectrum ^p{E) for a moder- 
ating medium with absorption; above Ep, but below maximal neutron energy E^""^ the 
spectrum was described by '^^^Pu fission spectrum [21, 22]. Here the neutron gas tem- 
perature for Maxwell distribution was given by (82), described in [35]. According to this 
approach [35], drawbacks of standard slowing-down theory for thermalization area may 
be formally reduced if a variable ^(x) = ^(1 — 2/z) is introduced instead of the average 
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logarithmic energy loss ^, which is almost independent of neutron energy (as is known, 
for environment consisting of nuclei with yl > 10 the statement ^ ~ 2/^ is true). Here 
z = En/kT, En is the neutron energy, T ~ environment temperature. Then within such 
a framework for neutron gas temperature in Maxwell spectrum of thermal neutrons the 
following expression may be used^: 



T 



l+T] 



ms 



(82) 



where Tq is the fuel environment temperature, Sa(A;To) is an absorption cross-section 
for energy /cTq, ry = 1.8 ~ dimensionless constant, (^) is averaged over the whole energy 
interval of Maxwell spectrum ^{z) at kT = 1 eV. 

Fermi neutron spectrum for a moderating medium with absorption (we considered 
carbon as a moderator and ^^^C/, '^^^U and ^^^Pu as absorbers) was set in the form 
[35, 54]: 



Fermi 



S 



exp 



Sg {E') dE' 

(OSt {E')E' 



(83) 



where S is the total volume neutron generation rate, (^) = Yl (^j^s) f^Si is the 

i 

average logarithmic decrement of energy loss, - macroscopic scattering cross-section 
of i-th nuclide, T,t = Yl '^s + ~ total macroscopic cross-section of the fissile material, 

i 

Tis = Y ~ total macroscopic scattering cross-section of the fissile material, - 

i 

macroscopic absorption cross-section, Ep - upper neutron energy for Fermi spectrum. 

The upper limit of neutron thermalization E^^ in our calculation was considered as 
a free parameter, setting the neutron fluxes of Maxwell and Fermi spectra at a common 
energy limit En^^ equal: 



Maxwell {Eum 



^Fermi {El 



lim ) 



(84) 



The high energy neutron spectrum part {E > Ep) was defined by fission spectrum 
[54, 65, 66] in our calculations. Therefore for the total volume neutron generation rate S, 
in the Fermi spectrum expression (83), the following expression may be written : 



S{r,T,t)= J P{r,E,T,t) 



u,{E) ■ $ (f, E, T, t) ■ a}{E, T) ■ Ni (f, T, t) 



dE, (85) 



where i?™"^ 



is the maximum energy of neutron fission spectrum (usually taken as EJ^""^ ~ 
10 MeV), Ep is the neutron energy, below which the moderating neutrons spectrum is 
described as Fermi spectrum (usually taken as Ep w 0.2 MeV); P (r, E, T, t) - probability 
of neutron not leaving the boundaries of fissile medium which depends on fissile material 
geometry and conditions at its border as well (e.g. presence of a refiector). 

The obtained calculation results show that cross-sections averaged over the spectrum 
may increase (fig. 12 for ^^^Pu and fig. 14 for '^^^U) as well as decrease (fig. 13 for "^^^U) 
with fissile medium temperature increase . As follows from the obtained results, the 
arbitrariness in selection of the limit energy for joining Maxwell and Fermi spectra does 
not significantly alter the character of these dependencies evolution. 

This can be justified by the fact that ^^^Pu resonance area starts from significantly 
lower energies than that of ^ss^y ^^j^^ with fuel temperature increase, the neutron gas 

very interesting expression revealing hidden connection between temperature of neutron gas and envi- 
ronment (fuel) temperature. 
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a) b) 

Figure 12: Temperature dependencies for fission cross-section (a) and radioactive capture cross- 
section (b) for "^^^Pu, averaged over the Maxwell spectrum, on the Maxwell and Fermi spectra 
joining energy and r] = 1.8 (see (82)). 
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Figure 13: Temperature dependencies for fission cross-section (a) and radioactive capture cross- 
section (b) for "^^^U, averaged over the Maxwell spectrum, on the Maxwell and Fermi spectra 
joining energy and rj = 1.8 (see (82)). 
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Figure 14: Temperature dependencies for fission cross-section (a) and radioactive capture cross- 
section (b) for 92^f/, averaged over the combined Maxwell and Fermi spectra depending on 
Maxwell and Fermi spectra joining energy and rj = 1.8 (see (82)). 
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temperature increases producing Maxwell's neutron distribution maximum shift to higher 
neutron energies. In other words, the neutron gas spectrum hardening, when more neu- 
trons fit into resonance area of '^^^Pu, is the cause of the averaged cross-sections growth. 

For [/this process in not as significant because its resonance area is located at higher 
energies. As a result, the 235[/ neutron gas spectrum hardening related to fuel temperature 
increase (in the considered interval) does not result in significant increase of a number 
of neutron fitting into the resonance area. Therefore according to the known expressions 
for 235^ determining the neutron reactions cross-sections behaviour depending on their 
energy En for non-resonance areas we observe dependency for the averaged cross-sections 



The data on the averaged fission and capture cross-sections of 238 presented at fig. 14 
on one hand show that the averaged fission cross-section for '^'^^U, due to a high fission 
threshold ~1 MeV (see fig. 14a), is almost insensitive to neutron spectrum hardening 
caused by the fuel temperature increase, and on the other hand - confirm the dependency 
of capture cross-section on temperature, because its resonance area is located as low as for 
239 Obviously in this case fuel enrichment with ^35^ essentially makes no diff'erence 
because the averaged cross-sections for '^^^U, as described above, behave in a standard 
way. 

And finally we performed a computer estimate of thermal source density dependence 
Qj, (r, $,T, t) (81) on temperature for different compositions of uranium-plutonium fissile 
medium with a constant neutron fiux density, presented at fig. 15. We used the dependen- 
cies presented above at fig. 12-14 for these calculations. Let us note that our preliminary 
calculations were made not taking into account the change in composition and density of 
the fissile uranium-plutonium medium that is a direct consequence of constant neutron 
flux assumption. 

The necessity of such assumption is caused by the following. It is obvious that for 
reasonable description of neutron source density qJp (r, ^,T,t) (81) temperature depen- 
dency a system of three equations must be solved, two of which correspond to neutron 
kinetics equation (flux and fluence) and system of equations for kinetics of nuclear density 
of parental and child nuclides (e.g. see [26, 29]), while the third corresponds to a heat 
transfer equation of (76) type. However some serious difficulties arise here, associated 
with the computational capabilities available . And here is why. 

One of the principal physical peculiarities of TWR is the fact [20] that fiuctuation 
residuals of plutonium (or 233[y Th-U cycle) over its critical concentration burnout 
for time comparable with reactor lifetime of neutron Tn{x,t) (not considering delayed 
neutrons) or, at least comparable with the reactor period^ T{x, t) (considering delayed 
neutrons), while new plutonium (or 233jy Th-U cycle) is formed in a few days (or a 
month) and not at once . This means [20], that numerical calculation must be made with 
temporal step around 10~^-10~^ in case of not taking into account the delayed neutrons 
and ~10~^-10'' otherwise. At first glance taking into account the delayed neutrons ac- 
cording to [20] really "saves the day", however it is not always true. If the heat transfer 
equation contains a significantly non-linear source, then in the case of a blow-up mode, 
the temperature at some conditions may grow extremely fast and in 10-20 steps (with 
time step 10~^-10~'^ s) reaches critical amplitude that may lead to (at least) a solution 
stability loss or (as maximum) blow-up bifurcation of the phase state, almost unnoticeable 
with rough time step. 

Here we should add that scientific literature has absolutely no theoretical or experimen- 
tal data on heat capacity Cp and heat conductivity bJ; of fissile material for temperatures 
over 1500K, which makes the further model calculations of the problem quiet problematic. 

According to these remarks, and considering the goal and format of this paper, we 
didn't aim at finding an exact solution of some specific system of three joined equations 

^The reactor period by definition equals to T(a;, t) — Tn{x, t)/ p{x, t), i.e. is a ratio of reactor neutron lifetime 
to reactivity. 
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Figure 15: Dependency of the thermal source density q;!p {r,^,T, Ni,t) [eV] on the fissile 
medium temperature (300-1400K) for several compositions of uranium-plutonium medium (1 - 
10% Pu; 2 - 5% Pu; 3 - 1% Pu) at constant neutron fiux density $ = 10^^ n/{crn? ■ s). 



described above. However, we found it important to illustrate - at qualitative level - 
the consequences of possible blow-up modes in case of non-linear heat source existence in 
heat transfer equation. As said above, we made estimate computer calculations of thermal 
source density qj, (r, <I>, T, t) (81) temperature dependence in 300-1400K range for some 
compositions of uranium-plutonium fissile medium at constant neutron flux (fig. 15). 

The obtained dependencies for source thermal flux density q^, (r, T, t) were success- 
fully approximated by a power function of temperature with an exponent of 4 (fig. 15). 
In other words, we obtained an equation for heat transfer with a significantly nonlinear 
heat source in the following form: 

qT{T)= const ■T^^+^\ (86) 

where 5 > 1 in case of non-linear thermal conductivity dependency on temperature [56- 
60]. The latter means that the solutions of the heat transfer equation (76) describe the 
so-called Kurdyumov blow-up modes [56-61], i.e. such dynamic modes when one of the 
modeled values (e.g. temperature) turns into infinity for a finite time interval. As noted 
before, in reality instead of reaching infinite values a relevant phase transition is observed 
( final phase of parabolic temperature growth), that requires a separate model and is a 
basis for an entirely new problem. 

Mathematical modeling of blow-up modes was performed mainly using Mathematica 
5.2-6.0, Maple 10, Matlab 7.0, utilizing multiprocessor calculations for effective applica- 
tion. For calculations a Runge-Kutta method of S-O*'' order and the numerical methods 
of lines [67] were applied. The numerical error estimate didn't exceed 0.01%. The coordi- 
nate and temporal step were variable and chosen by the program in order to fit the given 
error at every calculation step. 

Below we give the solutions for heat transfer equation (76) with nonlinear exponen- 
tial heat source (86) in uranium-plutonium fissile medium for boundary and initial pa- 
rameters corresponding to technical reactors. The calculations were done for a cube 
of fissile material with different spatial size, boundary and initial temperature values. 
Since the temperature dependencies of thermal source densities were obtained without 
account for changing composition and density of uranium-plutonium fissile medium, dif- 
ferent blow-up modes can take place (HS-mode, S-mode, LS-mode) depending on a ratio 
between the exponents of temperature dependencies of thermal conductivity and heat 
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Figure 16: Temperature dependency of heat capacity cp and thermal conductivity x of fissile 
material. Points - experimental values for heat capacity and thermal conductivity of "^^^U . 



source according to [56-61]. Therefore we considered cases for 1*^*, 2'^'^ and 4**^ temper- 
ature order sources. Here the power of the source also varied by varying the propor- 
tionahty factor in (86) {const = 1.00 J/(cm'^ • s ■ K) for 1*^* temperature order source; 
0.10J/(cm3 • s ■ K^), 0.15 J/{cm^ ■ s • K"^) and 1.00 J/{cm^ ■ s ■ K'^) for 2^^^^ temperature 
order source; 1.00 J /{cm? ■ s ■ K^) for 4*^ temperature order source). 

During calculations of the heat capacity Cp (fig. 16a) and thermal conductivity H 
(fig. 16b) of fissile medium dependency on temperature in 300-1400K range the specified 
parameters were given by analytic expressions, obtained by approximation of experimental 
data for^'^^C/ based on polynomial progression: 

Cp{T) ^ -7.206+0. 64r-0.0047r2+0.0000126T^+2.004-10"^r'^-1.60-10"^°r^-2.15-10"^^T^ 

(87) 

H(T) 21.575 + 0.015266ir. (88) 

And finally the thermal transfer equation (76) solution was obtained for constant 
thermal conductivity (27.5 W/{m ■ K)) and heat capacity (11.5 J/{K ■ mol)) values, 
presented in fig. 17a, and also the solutions of heat transfer equation considering their 
temperature dependencies (fig. 17b-d). 

Preliminary results point directly to a possibility of local melting of uranium-plutonium 
fissile medium, having melting temperature almost identical to that of that equals 

1400K (fig. 16a-d). Moreover, these regions of local melts are not the areas of so-called 
thermal peaks [68] and, probably, are the anomalous areas of uranium surface melt ob- 
served by Laptev and Ershler [69] that were also mentioned in [70]. More detailed analysis 
of probable temperature scenario associated with blow-up modes are discussed below. 

6 The blow-up regimes in neutron- multiplying me- 
dia and impulsive thermonuclear TWR. 

Earlier we noted the fact that due to coolant loss at nuclear reactors during Fukushima 
nuclear accident the fuel was melted or, in other words, temperature inside the active 
zone at some moment reached melting temperature of uranium-oxide fuel, i.e. ~3000°C. 

On the other hand, we already know that coolant loss may become a cause of non- 
linear heat source formation inside the nuclear fuel and therefore become a cause of 



29 



Figure 17: Thermal transfer equation (76) solution for 3D case (crystal sizes 0.001x0.001x0.001 
mm; initial and boundary temperatures equal to lOOK): a) The source is proportional to 4*^ 
order of temperature; const =1.00 J/(cm^ ■ s ■ K'^), heat capacity and thermal conductivity 
are constant and equal to 11.5 J/{K ■ mol) and 27.5 W/{m ■ K) respectively; b) the source is 
proportional to 4*"^ order of temperature; const =1.00 J /{cm? ■ s ■ K^)] c) The source is propor- 
tional to 2°^^ order of temperature; const =1.00 J /{cm? ■ s ■ K"^); d) the source is proportional 
to 2^^^^ order of temperature; const = 0.10 J/{cm^ ■ s ■ K"^). Note: in cases b)- d) heat capacity 
and thermal conductivity were determined by (87) and (88) respectively. 
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temperature and neutron flux blow-up mode occurrence. A natural question arises, if it 
is possible to use such blow-up mode ( of temperature and neutron flux) for initiation of 
certain controlled physical conditions under which the nuclear burning wave would reg- 
ularly "experience" a so-called controlled blow-up regime. It is quite difficult to answer 
this question definitely, because such fast process has some physical vaguenesses, any of 
which can become experimentally unsurmountable for such process control. 

Nevertheless such process is very elegant and beautiful from the physics point of view 
and therefore requires more detailed phenomenological description. Let us try to make it 
in short. 

As we can see from plots of capture and fission cross-sections evolution for ^^^Pu 
(fig. 12), the blow-up mode may develop actively at '^1000-2000K (depending on real 
value of Fermi and Maxwell spectra joining boundary), but at temperatures over 2500- 
3000K returns to almost initial cross-sections values. In case at that point we turn on 
some effective heat sink, the fuel may return to initial temperature state. However , while 
the blow-up mode develops, fast neutrons already penetrate to adjacent fuel areas, where 
new fissile material starts accumulating and so on (see cycles (1) and (2)). After some 
time in the this adjacent area the same blow-up mode will develop and everything starts 
over again. In other words, such hysteresis blow-up mode, closely time-conjugated to a 
heat takeoff procedure, will appear at background of stationary nuclear burning wave in 
a form of periodic impulse bursts. 

In order to demonstrate the marvelous power of such process, we investigated the heat 
transfer equation with non-linear exponential thermal source in uranium-plutonium fissile 
medium with boundary and initial parameters emulating the heat takeoff process. In 
other words, we investigated blow-up modes in fast Feoktistov-type uranium-plutonium 
reactor (1) where inside and outside boundary the temperature was deliberately fixed at 
6000K. 

This temperature is defined by an important question: "Is it possible in such conditions 
(6000K) emulating the time-conjugated heat takeoff (see fig. 12) to obtain a solution, i.e. 
spacio-temporal temperature distribution not in form of (5-function at some local spatial 
area, but as some kind of stationary and limited by amplitude solitary wave?" As shown 
below, such suggestion proved productive. 

Below we present some calculation characteristics and parameters. During these calcu- 
lations we used the following expression for dependency of thermal conductivity coefficient: 

H = 0.18 • lO"'* • T, 

which was obtained using Wiedemann-Franz law and data on electric conductivity of met- 
als at temperature 6000K [71]. Specific heat capacity at constant pressure was determined 
by value Cp « 6 cal/{mol ■ deg) according to Dulong and Petit law. 

Fissile uranium-plutonium medium was modeled as a cube with size 10.0 x 10. Ox 10.0 m 
(fig. 18) during calculations. Here for thermal source we used 2nd order temperature 
dependency (see (86)). 

And finally fig. 18a-d present a set of solutions of heat transfer equation (76) with 
nonlinear exponential heat source (86) in uranium-plutonium fissile medium with bound- 
ary and initial conditions emulating such process of heat takeoff that initial and boundary 
temperatures remain constant and equal to 6000K. 

It is important to note here, that solution set, presented at fig. 18, is q demonstrate 
quite clearly the temporal evolution of the solution to its "Stationary" state. This is 
achieved using the so-called "magnifying glass" approach when the solutions of the same 
problem are deliberately investigated at different timescales. E.g. fig. 18a shows the 
solution at time scale t E [0, 10~^ s], while fig. 18b describes spatial solution of the 
problem (temperature field) for t = 10~^ s. The fig. 18c-d presents the solution - the 
spatial temperature distribution - at t = 0.5 s and t = 50 s. 

As one can see, the solution (fig. 18d) is completely identical to the previous (fig. 18c), 
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Figure 18: Heat transfer equation solution for model reactor (source 2"*^ order temperature 
dependency, const = 4.19 J/{crrr' ■ s ■ K"^); Initial and boundary temperatures equal to 6000K; 
fissile medium is a cube 10x10x10 m. The presented results correspond to the following times 
of temperature field evolution: (a) (1-10)-10~'' s, (b) 10~^ s, (c) 0.5 s, (d) 50 s. 
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i.e. to the distribution established in the medium in 0.5 seconds, that allowed us to make 
a conclusion on temperature field stability, starting from some moment. It is interesting 
that the established temperature field creates the conditions enough for thermonuclear 
synthesis reaction, i.e. reaching lO'^K, and such temperature field existence time is not 
less than 50 s. With a corresponding nuclei concentration to enter thermonuclear synthesis 
reaction, these condition, according to a known Lawson criterion, are favorable for a stable 
thermonuclear burning. 

However, one should remember that the results of this chapter are only demonstrative, 
because their accuracy is very relative and requires careful investigations with application 
of necessary computational resources. However, qualitative peculiarities of these solutions 
should attract the researchers' attention to the nontrivial properties of blow-up modes, at 
least, with respect to the obvious problem of internal TWR safety violation. 

7 Conclusion 

Below we give short conclusions stimulated by significant problems, that can be formulated 
in the following form. 

1. TWR and problem of dpa-parameter in claddings materials. A possibility 
for surmounting a so-called problem of dpa-parameter based on conditions of nuclear 
burning wave existence in U-Pu and Th-U cycles. In other words, joined application 
of " differential" [29] and "integral" [20, 21, 26] conditions for nuclear burning wave 
existence allows to determine such nuclear burning wave mode, that has parameters 
(fluence/neutron flux, width and speed of nuclear burning wave) that completely 
satisfy the dpa-condition (73), that characterizes tolerable level of radioactive sta- 
bility of reactor materials, in particular of claddings materials. The latter means 
that at the present time the problem of dpa-parameter in claddings materials in the 
TWR-project is not an insurmountable technical problem and can be satisfactorily 
solved. 

Here we may add, that this algorithm of optimal nuclear burning wave mode se- 
lection predetermines a satisfactory solution of other technical problems mentioned 
in introduction. For example, the fuel rod length in a proposed TWR variant (see 
"ideal" case in Table 1) is predetermined by nuclear burning wave speed, which in 
a given case equals to 0.254 cm/day=85 cm/year, i.e. 20 years of TWR operation 
requires fuel rod length ~17 m. On the other hand, it is known [72], that for a 
twisted fuel rod form with two- or four-bladed symmetry, the tension emerging from 
the fuel rod surface cooling is 30% lower than that of a round rod with the same 
diameter, other conditions being equal. The same reduction effect applies to the 
hydraulic resistance value in comparison to round rod of the same diameter. 
Another problem associated with reactor materials swelling is also solved rather 
simply. Here we should note that if a ferritic-martensitic material is chosen as a 
cladding material (fig. 8 [48]), then the swelling effect at the end of operation will 
be only ^0.5% [48]. We could discuss other drawbacks noted in the introduction as 
well, but in our opinion , comparing to solution of a main problem predetermined 
by dpa-condition, the rest of the problems are not super-obstacles for contemporary 
level of nuclear engineering and can be solved in a traditional way. 

2. The consequences of anomalous 238fj ^^j^^j 239p^ cross-sections behavior 
with temperature increase. It is shown that the capture and fission cross-sections 
of 238f/ g^j^jj 239 manifest a monotonous growth in 1000-3000K range. Obviously 
such anomalous temperature dependency of '^^^JJ and ^^^Pu cross-sections cardinally 
change the neutron and thermal kinetics of nuclear reactors, TWR in particular. It 
becomes essential to know their influence on kinetics of heat transfer because it may 
become the cause of positive feedback with neutron kinetics, which may lead not 
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only to undesirable loss of solution stability (the nuclear burning wave), but to a 
trivial reactor runaway with a subsequent nontrivial disaster. 

3. Blow-up modes and problem of nuclear burning wave stability. One of 

the causes of possible temperature growth of the fuel is, for instance, a deliberate 
or spontaneous coolant loss analogous to Fukushima nuclear accident. As shown 
above, the coolant loss may become a cause of nonlinear thermal source formation 
in nuclear fuel, and therefore emerging of mode with temperature and neutron flux 
blow-up . In our opinion, the preliminary results of heat transfer equation with non- 
linear thermal source investigations point to an extremely important phenomenon of 
anomalous behavior of thermal and neutron flux blow-up modes. This result poses 
a natural nontrivial problem of fundamental nuclear burning wave stability and, 
correspondingly, of a physically reasonable application of Lyapunov method to this 
problem, which is a basis for motion stability theory, and thus a reliable basis for 
justification of the Lyapunov functional minimum existence. 

It is shown that some variants of solution stability loss are caused by anomalous 
nuclear fuel temperature evolution. They can be not only the cause of TWR internal 
safety loss, but through a bifurcation of states (and this is very important!) can lead 
to a new stable mode when nuclear burning wave would periodically "experience" 
the so-called controlled blow-up regime. At the same time, it is noted that such fast 
(blow-up regime) process has a number of physical vaguenesses, any of which may 
happen to be experimentally insurmountable for control of such process. 

4. On-line remote neutrino diagnostics of intra-reactor processes. Due to the 

fact that a high-power TWR or a nuclear fuel transmutation reactor are projects 
of single-load and fuel burn- up with a subsequent burial of the reactor apparatus, 
there is an obvious necessity for remote system of neutrino control of neutron-nuclear 
burning wave in normal operation mode and control of neutron kinetics in emergency 
situation. The details and peculiarities of calculation of spatial-temporal distribution 
of isotope composition of active zone in TWR are presented in [17, 73, 74] in detail 
within the inverse problem of neutrino diagnostics of intra-reactor processes. 
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